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Abstract 

We give the Green function, momentum distribution, two-particle correla- 
tion function, and structure factor for the bound state of N indistinguishable 
bosons with an attractive delta-function interaction in one dimension, and 
an argument showing that this boson "molecule" has no excited states other 
than dissociation into separated pieces. 
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I. THE BOSON MOLECULE 



Consider a system of iV identical bosons in one dimension, with an attractive delta- 
function interaction. There is a bound state for all N, having the form of a well-localized 
"molecule" JD-Q. This result is interesting, since it is the unique example of an exactly solved 
localized many-body system. A possible experimental realization of this system would be 
in a 7 Li vapor jl confined to an atomic trap which is so tight in two directions, that 
the system can be regarded as one-dimensional with respect to the third direction. Such 
traps are realistic prospects in the near future ||, and high-aspect ratio cigar-shaped traps 
approximating quasi-one-dimensional systems are already available experimentally ||. 

We seek the localized eigenfunctions \l/{xj} of the Hamiltonian 
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wherein m is the particle mass and W is a positive parameter representing the strength of 
interparticle attraction. 

The wavefunction and many derived quantities can be given exactly and in closed form; 
however, it is useful to first consider the mean-field version of this theory. The Gross- 
Pitaevskii (GP) approximation J/J assumes that ^{xj} can be written as a simple product 
of wavef unctions 

*(fo}) = 11^0*0 ( 2 ) 

i 

where the if)(x) are eigenfunctions for a single particle moving in the average potential of all 
the others 

= -l^-^r + V(x^(x). (3) 

where V(x) = — (N — l)W\ipo(x)\ 2 . Thus the ground-state wavefunction ijj is actually deter- 
mined by solving a nonlinear Schrodinger equation. We will show that this approximation 
generally captures the physics of the large N limit; this is in contrast to the case of repulsive 
bosons in one dimension, for which the interaction terms must be treated more carefully |§ . 
The solution to Eq.(3) is 



. . s (N-l)a 1 

^° (X) = V 4 cosh((iV-lW2) (4) 

where a = mW/h 2 ; the single particle energy is e = — a(N — l) 2 W/8. This has the impli- 
cations for the many-body problem that there can be a state that is localized in the sense 
that it is unlikely that any particle be found far from the center of mass, and provides an 
estimate for the particle density (with the center of mass fixed at the origin) n(x) = Nip(x) 2 ; 
the corresponding ground state energy is E = —aN(N — 1) 2 W/2A. 

It can be shown that there are no other bound states of Eq. (3) for the potential V(x) 
constructed from the ground state wavefunction. This suggests that we cannot construct ap- 
proximate localized excited states of the boson molecule by symmetrizing ip(xi) Y\^ =2 V'otX') 
(for some appropriate choice of ip). 
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The GP wavefunction is not translationally invariant, and thus violates momentum con- 
servation. Every translation of it is a solution to the same self-consistent equation with the 
same energy; then by forming the linear combination ^s{{ x i}) — I ^({^i — X})dX we can 
restore this symmetry. 

II. THE EXACT WAVEFUNCTION 

For the case N = 2, the wavefunction can be found by changing to relative coordinates, 
giving a problem equivalent to a single particle and a fixed potential well. The (unnormal- 
ized) ground state wavefunction is 

^2(^1,^2) = exp(--a\xi - x 2 \) (5) 

where a = mW/h 2 ; the corresponding energy is E 2 = —h 2 a 2 /Am = —mW 2 /Ah 2 . For 
x\ 7^ X2 the delta function vanishes, and the wavefunction is an eigenfunction of the kinetic 
energy operator; for a bound state this is the product of exponential functions with real 
argument. The potential generates the discontinuity in the derivative of $ at the crossing 
point x\ = x 2 . These two properties hold for all N, so that we can write down the general 
(unnormalized) wavefunction as || 

#jv(zi, • • • ,x N ) = exp(-(a/2)^2\xi - Xj\) (6) 

i<j 

For the sector x\ < x 2 < . . . < xn this can also be written in the form 

m N = exp((l - N)ax N /2 + (3 - N)ax N . 1 /2 + . . . + (N - l)ax 1 /2) (7) 

In view of the exchange symmetry of bosons, this completely describes the wavefunction. 
The ground state energy is proportional to the sums of the squares of the coefficients of X{ 
in the argument of the exponential; explicitly, 

E N = -(N 3 - N)aW/2A = — (iV 3 - N)mW 2 /2A% 2 (8) 

This agrees with the GP result only in leading order in N. 

III. THE SINGLE PARTICLE DENSITY 

The implications of the wavefunction are more clearly revealed when reduced to the 
probability that there is a particle at position x relative to the center of mass. For the case 
N = 2, there can only be a particle at x if the other is at —x, and then the probability is 
proportional to the square of the wavefunction: 

P 2 = 2aexp{-2a\x\) (9) 

This has been normalized so that the integral over all x gives N = 2. 
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For larger N we integrate over the unconstrained degrees of freedom; for example for 
N = 3 the probability of finding a particle at x > is proportional to 

/x poo 
^> 2 {-x - x 2 ,x 2 ,x)dx 2 + / ^> 2 {-x - x 3 ,x, x 3 )dx 3 (10) 
-x/2 Jx 

where the integrals describe the case that it is particle 3 or particle 2 that is at x (particle 
1 is necessarily in the region x < 0). The resulting normalized probability distribution is 

P 3 {x) = 3 • 2a[exp(— 3a|x|) — - exp(— 6a|x|)] = 3 ■ 2a[z — -z 2 ] (11) 

where in the second representation we have introduced the abbreviation z = exp(—Na\x\). 
In a similar way P/v can be constructed for other small N. The results for N < 7 suggest 
the general expression (given previously by Yoon and Negele 0) 

P„ W = a St" 1 ) (N + n _ 1)]{N _ n _ ijj eM-nNaW) (12) 

This has been normalized so that the integral over all x gives N. The density at x = is 
-Pjv(O) — N 2 (N — 1)/ {AN — 6). An interesting feature of the distributions is revealed when 
we represent them as a power series in \x\: the coefficients of |x| ra vanish when n is odd and 
less than 2N — 3. Thus the discontinuity in slope at x = that is present in P2 is replaced by 
a discontinuity in a much higher derivative for larger iV (see footnote 7). This property is a 
demonstrable consequence of the representation (12); in fact, combined with the assumption 
that Pjy is a polynomial in z of order {N — 1), this is sufficient to determine the form of P/v 
up to normalization. 

In the limit of large N, Eq. (12) reduces to 

Af_1 N 2 a 
P N {x) ss iV 2 a {-l) n+1 nexp{-nNa\x\) ss ^— sech 2 (iVaa;/2) (13) 

n=l ^ 

which has form similar to the GP result P/v(x) = a sech 2 ((iV — l)ax/2). Figure 1 

shows {\/N 2 a)Pjq{x) as a function of Nax, for various A. 

We may also characterize P/v(x) through its moments. Starting from (12) it can be 
shown that 

1 f°° 7V_1 2 7T 2 

<x2>= NiJ p " {x)ix = ^.^ww^i^m (14) 



jv-i N-i JV_1 12 77r 4 

<X >= £ £ a 4 AT 4 M 2 P 2 + ^ « 4 iV 4 M 4 ~ 15a 4 A^ 4 ^ 

The GP theory would give < x 2 >= 3q ,2(^_ 1 )2 (and similarly for < a; 4 >) - they agree 
at large N, but the exact result approaches the asymptotic limit from below, while GP 
approaches from above. 
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IV. GREEN FUNCTION AND MOMENTUM DISTRIBUTION 



The single-particle Green function is calculated from the wavefunction by construct- 
ing the product ^(x, x 2 , x 3 , ..xn)^*(0, x 2 , x 3 , . . . xn) and integrating out the spectator coor- 
dinates X2,x 3 , . . . ,xn- Explicit construction of the Gn for iV < 6 suggests the result 

N 

G n(x) = f£ E {2 + a\x\[N 2 - 1 - 4K(m)}}z K{m) (16) 

/iV 771=1 

where w = exp(— a\x\) and K(m) = [N 2 — 1 — (N + 1 — 2m) 2 ]/ A. With the choice = 1, 
Gn(x) is normalized so that Gn(0) = 1; a physically more meaningful normalization would 
be to choose J^ 00 G^{x)dx = 1, but we have not succeeded in finding the closed form 
representation for the corresponding g^. Since K(N + 1 — m) = K(m), the sum contains 
every term twice, except for m = (N + l)/2 when N is odd. 
In the large N limit, K(m) « (N — l)(2m — l)/2, and then 

Gjv « Nax/2 smh(Nax/2) = - / ^{x + y)Mv)dy (17) 

which is the Green function that one would construct from the translationally invariant GP 
wavefunction ^s- 

The momentum distribution is the Fourier transform of Gat 

1 N jy~ 2 isl / 2 

where k is the wavevector. Figure 2 gives a graph of NaG^ik) versus k/Na, for various N. 
As in the case of the repulsive boson system [|l], the momentum distribution is smooth: there 
is no Bose condensation. The width of the momentum distribution is set by the parameter 
a, so that the momentum distribution for the noninteracting case (a = 0) is attained in a 
natural way. In a highly localized wavefunction we would hardly expect to find long-range 
off-diagonal order! 



V. TWO-POINT CORRELATION FUNCTION 

The two-point correlation function is calculated from the squared wavefunction by choos- 
ing one of the particles to be Xi = and another to be at Xj = x (in all possible ways), and 
integrating out all other {x{\. We have chosen to normalize C^{x) so that the integral over 
all x gives unity. It proved harder to analyze, because its representation is more irregular: 
for example, 

C 6 (x) = — [-110w 5 + 64w 8 + 81w 9 + a\x\{A00w 5 + 256w 8 )} (19) 
30 

where again w = exp(—a\x\). Based on the explicit construction of CV for iV < 7, we find 

N-l 

C N (x) = m _ £{[(10K(m) -2N 2 ]K(m) + a\x\[N 2 - AK(m)]K(m) 2 }w K ^ (20) 
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where K(m) = mN — m 2 . Since K(N — m) — K(m), the sum contains every term twice, 
except for m = N/2 when N is even. The special value at x = (which is the probability 
of having two particles at the same place) is Cn(0) = (N + l)/6. Figure 3 gives a graph of 
Cn(x) /aN 2 for various N. In the limit of large N, we may approximate iV 2 ^> K(m) in the 
prefactors and then K(m) « Nm, with the result 

> , a r ._ , . cosh(iVaa;/2) 2 

Cjv i « TTri^a N 97 rrl (21) 

NK ' 4N X 1 sinh 3 (iVa|:r|/2) sinh 2 (iVaa;/2) 1 V ; 

The structure factor is given by the Fourier transform of Cn(x): 

_ a ™ 2r [10K(m)-2iV 2 ] [A 2 - 4if(m)][A-(m) 2 - k 2 ] 

CN[k) ~ N(N - 1) ^ ^ (m) 1 k 2 + K(m) 2 [k 2 + K(m) 2 } 2 1 (22) 



VI. EXCITED STATES 

In this section we will show that for any iV there is only one localized state with zero 
total momentum. There are other states of negative energy, but these can all be interpreted 
as uncorrelated smaller molecules. 

For indistinguishable bosons we need only specify the wavefunction for the sector x\ < 
X2 < • • • < %n> an d in the interior of this region (where no two coordinates are equal), the 
Hamiltonian operator reduces to the kinetic energy. The bound eigenstates of the kinetic 
energy have the form 

^({xih {h}) = exp(^ kiXi), (23) 

and the energy is determined to be E = — J2f=i h 2 k 2 /2m. To be normalizable it is necessary 
that J2iLi k% — (otherwise a uniform translation of all the particles could lead to indefinite 
growth of ). The wavefunction for the interacting system is a linear combination of 
degenerate functions of this form. At the sector boundaries (where Xj = Xj + \ for some 
j), the interaction gives rise to a discontinuity in the derivative of the wavefunction, which 
implies 

(d/dx j+1 - d/dxj^ttxi})^^. = -a^({xi}) lXj+1=Xj (24) 
This must hold for arbitrary choices of the remaining N — 2 variables; then in the sum of 



terms of the form (23) the only terms that are coupled by (^4j) are those for which the lists 
{ki} differ only at kj and kj_i - but then the two sum constraints on the {ki} ensure that 
there are only two such lists, which differ by having the values for kj and kj_i interchanged. 
Thus the general form of the wavefunction would appear to be constructable from a single 
list {ki} (which we will take to be ordered, so that ki > kj for all i < j) in the form 



N 



where the first sum is over all permutations of N objects, and P(i) is the i th member of the 
permuted list. The boundary condition (pi) implies relations among the coefficients Cp 
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(a - kpQj) + kp ij+1] )C P = -(a - k Q{j) + k Q[j+1) )C Q = -(a + k P(j) - k P{j+1) )C Q (26) 



(It should be noted that the permutations P and Q are related so that P(j) = Q{j + 1) and 
P(j + l) = Q(j)). 

The foregoing differs only slightly from the case of bosons with repulsive delta-potential 
interactions [JTJ] - the main difference is that the ki are real instead of purely imaginary. 
However, the remaining boundary conditions are different: for repulsive interactions, the 
wavefunction is delocalized and must be confined within a periodic box, while for attractive 
interactions, the space can be infinite if the wavefunction itself is localized. 

The relevant (and interesting) condition is that ^({xi}) should become exponentially 
small if the extremal particle at x^ is moved to large positive x, or if the particle at Xi 
is moved to large negative x. This requires that fcpm be positive and k P ^ be negative, 
for every term in fl25|) that has nonzero coefficient. The sum itself is over all permutations, 
which certainly will generate terms in which the k's have the wrong sign; but the coefficients 
will vanish (according to (|2T)|)) for all permutations such that in the list {kp^} it happens 
there are members k P ^ and fcp(j) such that i < j and kp(j) = k P ^ — a. The ground state 
wavefunction is a simple example of the working of this rule: the ordered list is of the form 
of a sequence (ki, k\ — a,k\ — 2a, k\ — 3a, . . .) and every possible permutation of it gives rise 
to a vanishing coefficient. 

The ground state is not the only possibility: we could have a list which is made of 
several distinct sequences, so that in one of its permutations it takes the form {k\,k\ — 
a, . . .k a — n a a, k b , k b — a, ..k b — n b a, k c , k c — a . . .). This gives some freedom in the choice 
for the starting elements k p and sublist lengths n p ; this generates the wavefunction for an 
excited state of the molecule. However, it does not represent a localized state. Suppose we 
simultaneously displace the particles {xi : i = 1,2, ..,n a + 1} a large distance to negative 
values. If Yh=i ^ < 0, this term will become exponentially large, giving a nonnormalizable 
wavefunction. However, if Y17=i K > there will be a different permutation of the list in 
which this sublist is at the end, again giving a nonnormalizable wavefunction when these 
particles are moved to large positive values. The remaining possibility is J27=t 1 k{ — - but 
now the wavefunction remains finite if it is separated into two widely separated pieces, one 
of which contains n a + 1 particles: this describes states in which the particles have formed 
several separate and uncorrelated molecules. 

In the foregoing argument we have assumed that the {ki} are real. Allowing them to be 
complex does not change the situation much: the argument above implies that the k, having 
a real part must be organizable into sublists as above, all having the same imaginary part; 
each sublist represents a separated molecule, now with finite momentum. From this point 
of view we can readily see that when two molecules collide, nothing happens - they pass 
through each other leaving each undisturbed 0. 
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FIGURES 
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FIG. 1. Particle density distribution. Here are displayed Pn{x)/o>N 2 as a function of Nax, 
for N = 2,3, 4, 5, 10, and oo. The cusp in P2 at y = is replaced by discontinuities in the higher 
derivatives for other N. 
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